Observation of the Doubly Cabibbo-Suppressed Decay K^K^n 
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We report the first observation of tlie doubly Cabibbo-suppressed decay Z)J — > K^-k~ using 
605 fb~^ of data collected with the Belle detector at the KEKB asymmetric-energy e'^e~ collider. 
The branching ratio with respect to its Cabibbo-favored counterpart B{Df — > K'^n~) /B{Df 
K'^K~n'^) is (0.229±0.028±0.012)%, where the first uncertainty is statistical and the second 
is systematic. We also report a significantly improved measurement of the doubly Cabibbo- 
suppressed decay K'^tt'^7v~, with a branching ratio B{D^ — > K'^tt'^tv~)/B{D'^ — > 
7i'-7r+7r+)=(0.569±0.018±0.014)%. 

PACS numbers: 13.25.Ft, 14.40.Lb, 11.30.Hv 



Cabibbo-suppressed (CS) and doubly Cabibbo- 
suppressed (DCS) decays play an important role in stud- 
ies of charmed hadron dynamics. CS decays of nearly all 
the charmed hadrons have been observed, while DCS de- 
cays have been observed for only the -D+ and D° mesons. 
The naive expectation for the DCS decay rate is of the or- 
der of tan"* 9c, where 9c is the Cabibbo mixing angle [H, 
or about 0.29% [1] relative to its Cabibbo-favored (CF) 
counterpart. Current measurements roughly support 
this expectation. It is natural to extend the searches for 
DCS decays to other charmed hadrons in order to further 
understand the decay dynamics of charmed hadrons and 
complete the picture. 

Furthermore, one expects that the branching ratio of 
D+ K+TT+TT- [i] is about 2 tan"* 9c relative to its CF 
counterpart since the phase space for _D+ if ~7r+7r+ 
is suppressed due to the two identical pions in the fi- 
nal state. This expectation is consistent with current 
experimental results Therefore, we also expect the 
branching ratio of K^K^tt^ is about 1/2 tan"* 9c 

relative to its CF counterpart. Lipkin ,5;j argues that 
SU(3) flavor symmetry \d\ implies 

B{D+ ^ K+K+t:-) B{D+ ^ K+n+i^-) . 

1 X — ; — ■ ■ — — = tan un, 

B{Dt K+K-TT+) B{D+ ^ K-TT+n+) 

(1) 

where differences in the phase space for CF and DCS 
decay modes cancel in the ratios. The above relation 
does not take into account possible SU(3) breaking effects 
that could arise due to resonant intermediate states in the 
three-body final states considered here [^. 

In this Letter, we report the first observation of the 
DCS decay K'^K^n" and its inclusive branching 

ratio relative to its CF counterpart, D'^ K'^ . 
We also report a new measurement of the inclusive de- 
cay rate K^tt^tt~ relative to its CF counter- 
part, —>■ K^ir^TT^. The current upper limit on 
B{D+ ^ K+K+TT-)/B{D+ K+K-TT+) is 0.78% at 
the 90% confidence level (C.L.) [7] and the world-average 
of the _D+ i^+TT+TT^ branching ratio is BiD'^ 
K+Ti+Ti-)/B{D+ /f-7r+7r+) = (0.68±0.08)% H. We 
also test the validity of prediction ^ . 



The data used in the analysis were recorded at the 
T(45') resonance with the Belle detector [1] at the e+e" 
asymmetric-energy collider KEKB . It corresponds to 
an integrated luminosity of 605 fb~^. 

and candidates are reconstructed using three 
charged tracks in the event. The initial event selection 
is similar to that in other Belle measurements. We re- 
quire that the charged tracks originate from the vicinity 
of the interaction point with impact parameters in the 
beam direction (z-axis) and perpendicular to it of less 
than 4 cm and 2 cm, respectively. All charged tracks 
are required to have at least two associated hits in the 
silicon vertex detector fio'l , both in the z and radial direc- 
tions, to assure good spatial resolution on the D mesons' 
decay vertices. The decay vertex is formed by fitting 
the three charged tracks to a common vertex and requir- 
ing a C.L. greater than 0.1%. Charged kaons and pions 
are identified requiring the ratio of particle identification 
likelihoods, Ck/{Ck + C-,^), constructed using informa- 
tion from the central drift chamber, time-of-flight coun- 
ters, and aerogel Cherenkov counters jTl|, to be larger 
or smaller than 0.6, respectively. In addition, we require 
that the scaled momentum of the charmed meson candi- 
date Xp = p*/\/0.25 • E'cM ~ be greater than 0.5 to 
suppress combinatorial background as well as D mesons 
produced in B meson decays. Here p* and i?cM are the 
charmed meson momentum and the total e^e~ collision 
energy, calculated in the center-of-mass frame, and M is 
the reconstructed invariant mass of the candidate. Fig- 
ure [1] shows the K-k-k and KKn [l3| invariant mass dis- 
tributions after the initial selections. The background 
levels are too high to observe DCS signals. 

We then apply further selection criteria, which are op- 
timized using real data samples since there are some dis- 
crepancies between the Monte Carlo (MC) simulation 
and the data in the relevant distributions. We use 10% 
of the data sample for optimization and the remaining 
90% for the measurement to avoid a possible bias when 
the same samples are used for both optimization and the 
measurement. Hence the final selection criteria are ob- 
tained in a blind manner. Assuming no signal in the DCS 
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FIG. 1: M{K-k-k) (top) and M{KKit) (bottom) distributions 
after tlie initial selections. The solid curve shows the CF 
decay channel and the dashed shows the DCS decay channel 
in each plot. 



decay channel, we maximize Afs/VT^B, where Ms is the 
CF signal yield which has similar properties to the DCS 
signal and Mb is the background yield from the sideband 
regions in the DCS sample. 

One of the selections related to the finite lifetime of 
charmed hadrons is the reduced (x^/d-o.f) for the hy- 
pothesis that the candidate tracks for the charmed meson 
decay products arise from the primary vertex. The pri- 
mary vertex is obtained as the most probable point of 
intersection of the meson's momentum vector and the 
e+e" interaction region. Due to the finite lifetime of D'^ 
and mesons their daughter tracks are not likely to be 
compatible with the primary vertex. The second require- 
ment uses the angle between the charmed meson momen- 
tum vector, as reconstructed from the daughter tracks, 
and the vector joining its production and decay vertices. 
In an ideal case without resolution the two vectors would 
be parallel for the signal. The reduced is required to 
be greater than 25 for D+ and 5 for candidates and 
the angle is required to be less than 1° for and 2° for 
Df candidates. Tighter requirements on charged kaon 
identification (>0.9) and Xp (>0.7) are also chosen for 
the final selection, which improves the signal sensitivity. 
After the additional and tighter selection requirements 
described above, 9.57% of D+ and 10.71% of D+ CF sig- 
nal, and 0.06% of D+ and 0.24% of D+ DCS background 
events are retained. In order to minimize systematic ef- 
fects we choose the same selection criteria for both DCS 
and CF decay channels. The Ktttt and KKir invariant 
mass distributions after the final selections are shown in 
Figs, m and [3] together with signal and background pa- 
rameterizations. A clear signal is observed in both DCS 
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FIG. 2: Distributions of A/(7i'~7r+7r+) (top) and 
M(JC+7r+7r") (bottom). The K~-!t+tv+ distribution is 
shown on a semi-logarithmic scale to make the small back- 
ground visible. Points with error bars show the data and 
histograms show the results of the fits described in the 
text. Signal, Df — » K'^K'n^ background, and random 
combinatorial background components are also shown. 

The KiTTT and KKt: invariant mass distributions are 
fitted using the binned maximum likelihood method. In 
all cases the signal probability density function (PDF) is 
parameterized using two Gaussians with the same central 
value. Due to K/tt misidentification the following reflec- 
tion backgrounds appear in the mass distributions. In 
D+ -> K~TT+TT+ (CF) and K+TT+7r~ (DCS) decays there 
is a contribution from misidentified K^K~tt^ 
decay; in Z)+ K^K^n^ (CF) decay there is a con- 



tribution from misidentified 



and in 



£)+ if+i^+TT" (DCS) decay there is a contribution 
from misidentified K^tt^tt^ decay. The PDFs 

for the reflection backgrounds are determined from real 
data by assigning the nominal pion (kaon) mass to a kaon 
(pion) track. The magnitude of each of the reflection 
background contributions is a free parameter in the fit. 
For the DCS Df channel, the D+ i4:+7r+7r" contribu- 
tion is not incorporated in the fit since it is not signifi- 
cant, but its effect is included as a systematic uncertainty 
due to fitting listed in Table HIl The contribution 
{D*+ D°n+ with D" ^ K+K-) in the CF D+ chan- 
nel is also incorporated in the CF fit as an indepen- 
dent Gaussian component. A linear function is used for 
the random combinatorial background for all channels. 
All signal and background parameters for the CF chan- 
nels are fioated. For the DCS channels the mass, width, 
and ratio of the two signal Gaussians are fixed to the val- 
ues obtained from the fits to distributions of CF decays. 
Signal and background yields are left free in the fit. From 
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FIG. 3: Distributions of M{K+K-n+) (top) and 
M{K'^ K^TT^) (bottom). Points with error bars show the 
data and histograms show the resuhs of the fits described 
in the text. Signal, D*"*" background (Z)*"*" D^tv^ where 
D" — > K^K~), D"*" — > K~-k^tt'^ background, and random 
combinatorial background components are also shown. 



the results of the fits, shown in Figs. [5] and [31 we extract 
the signal yield for each channel, listed together with the 
corresponding branching ratios in Table [H 

The statistical significance of the Df K^K^n^ 
signal is calculated using —2\n{Ch/Cs+b) where Cb and 
Cs+b are the likelihood values of the fit, without and 
with the signal PDF included, respectively. We find 
— 21n(£b/£s+b)=83.2 with one degree of freedom used 
to describe the DCS signal yield; we obtain a statistical 
significance corresponding to 9.1 standard deviations. 

In addition to the backgrounds mentioned above there 
is also the possibility of double misidentification leading 
to contributions from CF events to the DCS sample. MC 
simulation shows that such a contribution is flat in the 
invariant mass distribution and is hence included in the 
combinatorial background description. 

The final states in this study have resonant substruc- 
ture that can affect the reconstruction efficiency. The 
resonances are relatively well known for the decay modes 
other than Z?+ K^K^n^ . We used a coherent 
mixture of resonant contributions according to [l4| to 
generate — > if~7r+7r+ decays and calculate the re- 
construction efficiency. For the £>+ K^tt^tt^ and 
Df K^K^t:^ decays we used an incoherent mix- 
ture of intermediate states Q. Subsequently we var- 
ied the contributions of individual intermediate states 
in a correlated manner, within the uncertainties of the 
measured branching fractions. The efficiency calculated 
from the modified MC sample differs from the original 
one by 1.5% and 2.0% for the D+ K+tt+tt- and 
Df K^K~n~^ decays, respectively, and the differ- 



ence was included in the systematic uncertainty of the 
result. Df K^K~^tt~ decays were generated accord- 
ing to phase space. For comparison, signal events were 
generated assuming either K*°{892)K+ or K*'^{USO)K+ 
intermediate states. The largest relative difference in the 
efficiency (2.4%) was included as a part of the systematic 
uncertainty. Ratios of reconstruction efficiencies for DCS 
and CF decays are found to be 1.042±0.008±0.016 and 
0.963±0.010±0.030 for D+ and Df decays, respectively, 
where the first uncertainty is due to the finite MC simu- 
lation statistics and the second is the uncertainty in the 
resonant structure of the final states. 

With the efficiencies estimated above, we measure the 
inclusive branching ratios of DCS decays relative to their 
CF counterparts summarized in Table HI The product 
of the branching ratios for the two DCS decay modes is 



found to be 



BiDt^K+K-7v+) ^ B(D+^K- 



1.57± 



0.21) X tan 9c, where the error is the total uncertainty. 



TABLE I: Measured branching ratios. Broi. is the branch- 
ing ratio relative to — > K~it'^-k^ for the D"*" modes and 
K~-K^ for the Dt modes. The uncertainties in the 
branching ratios are statistical and systematic. 



Decay Mode 



D+ 



signal 



S,el.(%) 



K- 



TT 



2637.7±84.4 0.569±0.018±0.014 
482702±727 100 



Di 



K 
K 



K 



TT 281.4±33.8 0.229±0.028±0.012 
7r+ 118127±452 100 



Several sources of systematic uncertainty cancel in the 
branching ratio calculation due to the similar kinematics 
of CF and DCS decays (for example, uncertainties in the 
tracking efficiencies and particle identification, since the 
momenta of the final state tracks are almost identical). 
The stability of the branching ratios against the varia- 
tion of the selection criteria was studied and we observed 
no changes greater than the expected statistical fluctua- 
tions. The systematic uncertainties due to the variation 
of the fit parameters are 1.9% for D"*" and 4.2% for 
branching ratios measurements, respectively. Table HI] 
summarizes the systematic uncertainties in the measure- 
ments of the branching ratios. 



TABLE II: Relative systematic uncertainties in percent, 

j^+j are systematic uncertainties 



where o"b^,^| {d+) ^^'^ '^b 



for the branching ratio of Z)+ and Df DCS decays relative to 
their CF counterparts. 



Source CTg 






Fitting 


1.9 


4.2 


MC Statistics 


0.8 


1.0 


Reconstruction Efficiency 


1.5 


3.1 


Total 


2.5 


5.3 
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Using the world-average values B{D~^ K 7r+7r+) = 
(9.22 ± 0.21)% and B{D+ K+K~Tr+) = (5.50 ± 
0.28)% 0, we obtain the absolute branching fraction for 
each DCS decay channel. Table Hill shows the results and 
the comparison to previous results. 



TABLE III: Absolute branching fraction for each decay mode 
and comparisons with previous measurements. The first un- 
certainties shown in the second column are the total uncer- 
tainties of our results and the second are the uncertainties in 
the average CF branching fractions used for normalization 
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Branching Fraction 



Belle 



B{Dt ^ K+K+Tv- 



World- Average [3] 
1 — 



(5.2±0.2±0.1)xl0"'' (6.2±0.7)xl0 
(1.3±0.2±0.1) X 10"" (2.9±1.1) X 10"^ 



To conclude, using 605 fb^^ of data collected with the 
Belle detector at the KEKB asymmetric-energy e+e" 
collider we have observed for the first time the decay 
K^K^TT^ with a statistical significance of 9.1 
standard deviations. This is the first DCS decay mode 
of the Z?^ meson. The branching ratio with respect 
to the CF decay is (0.229±0.028±0.012)%, where the 
first uncertainty is statistical and the second is sys- 
tematic. We have also determined the Z)"*" DCS de- 
cay branching ratio, B{D^ — > K^h^tt^)/B{D^ 
fs:-7r+7r+)=(0.569±0.018±0.014)%, where the first un- 
certainty is statistical and the second is systematic, with 
a significantly improved precision compared to the cur- 
rent world-average Q . We find the product of the two rel- 

B(D+^K-^K+Tr-) B(_D + ^_R-+7i-+7r~) 



ative branching ratios. 



B{Dt^K+K-T,+ ) B(D + ^K--K+-n+) 



to be (1.57±0.21)x tan**6'c. This is consistent with 
SU(3) flavor symmetry within three standard deviations; 
note that the effect of (different) resonant interniediate 
states is not taken into account in the prediction [3] . An 
amplitude analysis on a larger data sample may allow a 
more precise test of SU(3) flavor symmetry to be per- 
formed. 

We thank the KEKB group for excellent operation of 
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